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Abstract 

Microcanonical calculations are no more difficult to implement than canon- 
ical calculations in the Lattice Gas Model. We report calculations for a few 
observables where we compare microcanonical model results with canonical 
model results. 
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I. INTRODUCTION 



The Lattice Gas Model (LGM) is frequently used to compute observables in heavy ion 
collisions. The applications are numerous; see |T],|2[] for references. The calculations are 
normally done for a fixed temperature so that we will also call the the usual model canonical 
lattice gas model CLGM. The temperature is fixed so there will be fluctuations in energy. 
It is often argued that it might be more appropriate to keep energy fixed. In the following 
we set up a scheme for doing calculations with fixed energy. We will call this MLGM. 



II. CALCULATIONAL PROCEDURE 

We assume the reader is familiar with the usual LGM model for nuclear disassembly. 

In CLGM, N neutrons and Z protons are put in N s lattice sites using a Metropolis 
algorithm. Because of bonds between nearest neighbours [e np = — 5.33MeV, e nn = e pp = 0, 

and coulomb interaction between protons, there is a potential energy which we denote by 
E pot . In Metropolis method, switch is attempted between occupied sites and unoccupied sites 
and also between occupied neutrons and protons. If, in the switch, the energy goes down the 
move is accepted. If the energy goes up, the move is accepted but only with a probability 
exp(— AE/T). After many such switches an event is chosen. Once an event is picked, 
momenta are assigned from Monte-Carlo sampling of a Maxwell-Boltzmann distribution at 
temperature T. We denote the kinetic energy by E kin . The total energy of the system is 
E to t = Ep t + Ekm which will fluctuate from event to event. 

The change to MLGM can be made in the following way. We start from a given config- 
uration, hence a given E pot . The total energy which will be kept fixed is E to t, so the kinetic 
energy is E kin = E tot — E pot . The available phase space of A = N + Z nucleons having this 
kinetic energy is known analytically: 

/ S(E kin - f R 2 /2m)n^ = |^(2^) 3A/2 ~ 1/2 (2-1) 
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We will call the value of the integral Qi(Eki n )- We now attempt a switch in the configuration 
space. The potential energy will change to E ' t . To conserve energy the kinetic energy of 
the system should be fixed at E to t — E' pot = E' kin . If Qi(E kin )/fli(Eki n ) > 1, the move is 
accepted. If it is less, the probablity is given by the ratio. Since all configurations have 
identical weights, this satisfies the principle of detailed balance. After many such switches 
an event is accepted. We now have to assign momenta to the nucleons so that the total 
kinetic energy of the A nucleons add up to correct kinetic energy which we denote by Eki n . 
The correct way to do this such that the sole criterion is phase space is the following. Choose 
a sphere of radius P. Do a Monte-Carlo sampling on A particles for uniform distribution in 
this sphere. This means fixing p, 9 P , <p p for each particle from p = P(xi) 1 ^ 3 , cos9 p = 1 — 2x 2 
and (p p = 2ttx3 where X\,X2 and £3 are random numbers. Finally normalise P so that the 
total energy equals E kin . We are now ready to calculate all relevant quantities including 
cluster decomposition. 

In microcanonical simulation observables can be calculated without having to invoke a 
temperature. But it is useful to extract a temperature in the model. The "temperature" will 
vary from event to event. The event temperature is taken from \{A — 1)T rs \AT = E kin . 
The ensemble average gives the average temperature for the given microcanonical total 
energy. This is obviously attractive from an experimental point of view but we can also 
justify it from more basic principles. Let us write 



Although we are talking of one system only which has both kinetic and potential energy, 
formally, the right hand side is the same as two systems in "thermal" contact whose total 
energy is fixed but each one's individual energy can vary. This is very standard statistical 
mechanics [|J. For large systems, the sum is dominated by the maximum in the product 



nAEki^iEtot - E km ) which is obtained when dln ^^ = ^^i) which then defines 

unjoin O-Cypot 



the inverse of the average temperature. This of course leads again to the same identifica- 
tion as above and is consistent with eq.(2.1). The sharpness of the maximum in the sum 




(2.2) 
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depends on the size of the system. For a small system we will then expect the fluctuation 
in temperature < T 2 > — < T > 2 to be larger for the same value of < T >. 

Finally although the starting point of Metropolis algorithm can be fairly arbitrary, it is 
helpful computationally to start from close to equilibrium. The method of generatiing the 
starting point followed the scheme given in ||] which was modified to take into account that 
there are two kinds of bonds and the coulomb force between protons. 

Below we consider a few applications. Many more can be made. 

III. THE CALORIC CURVE 

As the first application, we show in Fig.l, the caloric curve, E*/A against temperature 
T, for 84 Kr (an intermediate mass disintegrating system) and for 197 Au (a heavy system). 
For the microcanonical model, < T > is used for T and for the canonical one, < E* /A > is 
used for E*/A. Even though the fluctuation < T 2 > — < T > 2 is more than four times in 
the case of Kr (6.59 MeV 2 at < T >= 5.0 MeV as compared to 1.1 MeV 2 for Au at the same 
average temperature) there is almost no difference in the E* j A vs. T curves. We should 
point out, all calculations use two kinds of bonds and include the coulomb interaction. For 
84 Kr we use a cubic lattice of size 6 3 ; for 197 Au, we used 8 3 . 

IV. IMF EMISSION PROBABILITY 

One might expect that while average quantities will be nearly the same in both the 
models, there would be larger differences in fluctuations of observables. With this in view, 
we have investigated Intermediate Mass Fragment {IMF, Z between 3 and 20) emissions. We 
show in Fig. 2, the plots of o 2 =< n 2 IMF > — < n IMF > 2 (upper panel) and a 2 / < n IMF > 
(lower panel) at different temperatures. Not a great deal of difference is found between 
CLGM and MLGM calculations for intermediate or heavy masses. Since the fluctuation in 
temperature for the intermediate case is much larger (more than a factor of four) this calls 
for an explanation. We think this is the reason. Referring to eq. (2.2) we see that when the 
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temperaure is low, E pot is high. That means in the particular configuration the number of 
attractive bonds is less. However, the probability of an attractive bond being able to bind 
two nucleons is much higher since the temperature is low (the energy of relative motion 
p 2 /2fi has a lesser chance of exceeding the bond energy — e np , see Similar arguments 
hold when the temperature is high. This means E pot is low, so that there are many more 
attractive bonds. However, higher temperature will be able to break these bonds more easily. 
The two effects seem to cancel each other quite efficiently. 

In Fig. 3 we compare P(tiimf), the probability of emitting n IMF's, in the two models 
for the case of 197 Au for four different temperatures. Again, the results are quite close. 

We verified that in Fig. 3, for temperatures 4.5 MeV and 5 MeV both the microcanon- 
ical and canonical probabilities are fitted quite well by a binomial distribution: P(n) = 
n\(m-n)\ P n (^ ~ p) m ~ n where m and p are obtained from < n >= mp and a 2 / < n >= 1 — p. 
This is a topic that has been discussed thoroughly in recent times [0. This parametrisation 
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does not work for temperatures 3.5 MeV and 4 MeV where -^—^ > 1. Here our calculated 
points are fitted quite well by a negative binomial distribution: P(n) = ^^1^ p n (l — p) N . 
The increase of a 2 at 4 MeV is not unnatural in LGM. This is happening because at this 
temperature the system is crossing the co-existence curve || with the accompaniment of a 
maximum in the fluctuation. 

V. SUMMARY 

To summarise, microcanonical calculations in LGM are no harder to do than canonical 
calculations. The lowest mass number that we used was 84. Down to this size at least there 
are no serious departures from canonical results. 
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FIGURES 



FIG. 1. Caloric curves for 84 Kr and 197 Au systems in the two models. 

FIG. 2. a 2 =< n 2 IMF > — < ujmf > 2 (upper panel) and a 2 / < ujmf > (lower panel) at 
different temperatures, as obtained in the two models. The left panels are for m Kr and the right 
ones are for 197 Au. 



FIG. 3. P(n IMF ) distributions for 197 Au at T = 3.5,4.0,4.5 and 5.0MeF in the two models. 
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